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Abstract 

The paper [l| by Crouseilles, Einkemmer, and Faou used an incorrect Poisson 
bracket for the Vlasov-Maxwell equations. If the correct Poisson bracket is 
used, the solution of one of the subsystems cannot be computed exactly 
in general. As a result, one cannot construct a symplectic scheme for the 
Vlasov-Maxwell equations using the splitting Hamiltonian method proposed 
in Ref. [ij. 
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In a recent paper [l| by Crouseilles, Einkemmer, and Faou, a new sym¬ 
plectic splitting method for the Vlasov-Maxwell equations is proposed. In 
comparison with previous splitting methods, the exciting new feature of the 
proposed method is that it is designed to preserve the symplectic structure 
of the Vlasov-Maxwell system, and thus enjoys the benehts of symplectic in¬ 
tegration, such as the global bound on energy error and long-term accuracy 
and fidelity. 

Crouseilles, Einkemmer, and Faou developed an innovative technique to 
achieve this goal. A non-canonical Poisson bracket for the Vlasov-Maxwell 
system as an infinite-dimensional Hamiltonian system is employed. The sys¬ 
tem is split into three subsystems by splitting the Hamiltonian functional 
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into three parts. It turns out that the solution for each subsystem can be 
computed exactly and therefore preserves exactly the symplectic structure 
corresponding to the Poisson bracket. As a consequence, the combined al¬ 
gorithm according to the splitting scheme preserves the symplectic structure 
as well. In addition, higher order methods can be constructed using various 
familiar composition methods. 

The Poisson bracket adopted by Ref. [l| is the bracket discovered by Mor¬ 
rison in 1980 j^. 
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for functionals E and G of E, B, and /. Here {h, is the canonical Poisson 
bracket in the (x, v) space. The Hamiltonian for the system is 


J {E^ + B‘^)dx. ( 2 ) 

This Hamiltonian can be split into three parts [H as follows. 


H = Hf 

+ He + Hb, 
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The scheme developed in Ref. [ij is based on the observation that solutions 
of the subsystems corresponding to Hf, He, and Hb can all be computed 
exactly. 

Unfortunately, this Poisson bracket ([I]) is known to be incorrect, because 
it does not satisfy the Jacobi identity. This error had been discovered and 
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corrected j^, 1^ shortly after its publication j^. The correct Poisson bracket 

0,8 is 
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Following Ref. 0]) will call this bracket the Morrison-Marsden-Weinstein 
(MMW) bracket. Integrating the third term on the right-hand side of Eq. ([7]) 
and considering the fact that 
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we can recast the MMW bracket as 
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In Ref. 0, the MMW bracket in the form of Eq. ([7]) is used, and in Refs. 0, 
0, the equivalent form of Eq. (|9]) is used. We emphasize that in order for 
the MMW bracket in the form of Eq. (JTj) or Eq. ([9]) to satisfy the Jacobi 
identity, some constraints in terms of B and/or E are necessary. Marsden 
and Weinstein 0 restricted the solution space to be Mv = {{f,E,B) \ 
V-R = 0, V-F = f fdv}. Morrison 0 pointed out that it is only necessary 
to require V • R = 0 for the MMW bracket in the form of Eq. ([7]) or Eq. (|9]) 
to satisfy the Jacobi identity. Of course, a solution of the Vlasov-Maxwell 
equations is always in Mn, if it is initially in Mv. 


3 








If we use this correct Poisson bracket instead, in the form of either Eq. ([7]) 
or Eq. ([2]), and apply the same splitting scheme as in Eq. (E]), it is very disap¬ 
pointing to hnd out that the solution for the subsystem corresponding to Hf 
can not be computed exactly, whereas the solutions for the subsystems cor¬ 
responding to He and Hb can. Specihcally, the subsystem associated with 
Hf is 
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This subsystem is more complicated than its counterpart obtained using the 
incorrect Poisson bracket ([H). Unless the magnetic held B is uniform in space 
or vanishes, its solution cannot be computed exactly. As a result, for systems 
with general magnetic held, one cannot construct a symplectic scheme for the 
Vlasov-Maxwell equations using the splitting Hamiltonian method proposed 
in Ref. [l|. 

If a symplectic integration method to a desired order for Eqs. flTUll - flTTD 
can be found, then we can apply this splitting method to obtain a symplectic 
scheme. However, such a symplectic method for Eqs. ffT0|) - ffT^ is not available 
yet. Further investigation is needed. 

As a hnal note, it is necessary to mention the following bracket proposed 
by Chandre et ah to remove the V ■ R = 0 constraint for the MMW 
bracket for the Vlasov-Maxwell equations. 
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Here, B — VA ■ R is the projection of R that “removes” the non¬ 
divergence-free part of R. It is straightforward to verify that the Jacobi 
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identity is unconditionally satisfied jsj. For this bracket given by Eq. flT^ . 
the splitting Hamiltonian method proposed in Ref. [l| generates identical 
subsystems as for the MMW bracket given by Eq. ([7]) or Eq. ([9]), when the 
V-R = 0 constraint is satisfied initially. Therefore, the splitting Hamiltonian 
method proposed in Ref. [ij is not valid for the bracket given by Eq. flT^ as 
well. 
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